Designing Symmetric Gradient Honeycomb Structures with Carbon-Coated Iron-Based Composites for High-Efficiency Microwave Absorption

Highlights MIL-88C (Fe) with varying aspect ratios as a precursor was synthesized by regulating oil bath conditions, followed by one-step thermal decomposition to obtain carbon-coated iron-based composites. High-efficiency microwave absorption properties were achieved with RLmin value of -67.4 dB (2.13 mm) and wide effective absorption bandwidth (EAB) of 5.52 GHz (1.90 mm) under the low filler loading. A symmetric gradient honeycomb structure was constructed utilizing the high-frequency structure simulator, achieving an EAB of 14.6 GHz and a RLmin of -59.0 dB. Supplementary Information The online version contains supplementary material available at 10.1007/s40820-024-01435-z.


Introduction
In the contemporary era, electromagnetic (EM) technology is undergoing continuous development and finds widespread applications in both military and civilian domains, including military radar, portable electronics, wireless communication, and more [1][2][3].While EM technology provides substantial convenience, it also gives rise to electromagnetic interference and radiation, posing risks to people's wellness [4][5][6], and impacting the regular operation of electronic devices.To address the aforementioned challenges, the development of materials that can absorb electromagnetic waves, particularly microwaves, is an imperative.Microwave absorbing materials (MAM) can be categorized into three groups based on distinct absorption mechanisms: materials with dielectric loss, magnetic loss, and composite loss [7].Dielectric loss materials primarily encompass carbonbased materials [8][9][10], MXene [11], non-magnetic metal compounds [12], and nonmetallic compounds [13].Notably, carbon-based materials have garnered significant attention owing to their distinctive properties, including lightweight characteristics, suitability for practical engineering applications, and favorable electrical properties.Magnetic loss materials predominantly consist of magnetic metals [14] and their compounds [15].Adequate impedance matching is a fundamental requirement for optimal microwave absorbing material (MAM) performance.Individual dielectric loss materials or magnetic loss materials alone cannot achieve effective impedance matching.Combining both types can overcome the aforementioned limitations [16].
Metal-organic frameworks (MOF) are coordination polymers comprising metal ions coordinated with organic ligands, forming intricate structures with diverse functionalities and applications.Composites comprising magnetic metals or their compounds coated with carbon can be synthesized by controlling pyrolysis conditions.These composites effectively combine dielectric losses from carbon materials and magnetic losses from either magnetic metals or their compounds, thereby enhancing impedance matching and MA performance [17].Currently, several studies on iron-based MOF derivatives with abundant metal element reserves and low costs have been conducted.Among the iron-based MOF, the MIL-88 series comprises four subtypes: 88A (fumaric acid), 88B (terephthalic acid, 1,4-BDC), 88C (2,6-naphthalenedicarboxylic, 2,6-NDC), and 88D (4,4′-biphenyldicarboxylic, 4-4′-BPDC), distinguished by different organic ligands [18,19].Among them, MIL-88C is characterized by a short synthesis cycle and a large AR.According to percolation theory, the larger the AR of one-dimensional (1D) materials, the lower the percolation threshold value.This scenario is conducive to achieving optimal MA performance at low filler levels [20].Previous studies have demonstrated that different ligands can influence the EM parameters of the corresponding derivatives, thereby impacting MA performance [21].Wu et al. [22] utilized MIL-88A as a precursor to synthesize the absorber Fe 3 O 4 /Fe 3 C@C through pyrolysis at 600 °C.The RL min was -52.9 dB at 11.68 GHz (3.07 mm) and a higher filler loading of 40 wt%, achieving a EAB of 4.56 GHz at 3.5 mm.Xiong et al. [23] synthesized Fe 3 O 4 @nanoporous carbon (NPC) composites by calcining Fe-MIL-88B.With a high filling amount of 30 wt%, Fe 3 O 4 @NPC exhibited a RL min of -65.5 dB at 9.8 GHz, while exhibiting an undesirable EAB of 4.5 GHz (3 mm).All of the aforementioned absorbers exhibited high filler amounts, narrow EAB, and thick matching thicknesses.Additionally, current research has focused on investigating the impact of pyrolysis conditions on MA properties.The influence of varying preparation conditions of the precursors on the MA properties of derivatives has been scarcely investigated.
The potential for practical engineering applications of single MAM is somewhat limited by inherent constraints.Numerous materials with unique structures exhibit MA properties surpassing those of natural materials, especially EAB [24][25][26].The honeycomb structure has gained prominence in current research due to its lightweight and highstrength characteristics [27,28].Impedance matching and attenuation capabilities are pivotal considerations in determining MA properties.Conventional through-hole honeycomb structures possess certain MA properties but do not independently isolate the two aforementioned factors.The entire honeycomb through-hole is responsible for both the entry and absorption of microwaves.Designing a structure that optimizes these two factors separately can more effectively leverage structural advantages and exhibit enhanced MA properties.
Herein, we controllably synthesized MIL-88C (Fe) as a precursor with varying AR using an oil bath under diverse conditions, followed by one-step thermal decomposition to obtain the corresponding carbon-coated iron-based composites.The influence of MIL-88C (Fe) preparation conditions 1 3 on the MA properties of the derivatives was investigated, including M/O, oil bath temperature, and oil bath time.Consequently, the phases, graphitization degree, and AR of the derivatives were influenced to attain low filler loading, optimize impedance matching, and enhance MA performance.As anticipated, the composites derived from MIL-88C (Fe) fabricated under a diverse range of conditions consistently displayed outstanding MA properties, including low loading, strong absorption, broad EAB, and thin thickness.We also designed a SGHS that concentrated on optimizing impedance matching in the upper half to facilitate greater microwave entry and on enhancing attenuation capability in the lower half to enable more efficient microwave absorption.Excellent MA performance was achieved through HFSS EM simulation.

Synthesis of MIL-88C (Fe)
MIL-88C (Fe) was synthesized using an oil bath method, a modified version of a previously reported fabrication procedure [29].Initially, 8 mL of DMF was stirred magnetically for 10 min with 75 mg (0.18 mmol) Fe(NO 3 ) 3 •9H 2 O and 4 mg (0.018 mmol) 2,6-NDC (M/O = 10:1).Subsequently, the mixture was subjected to heating in an oil bath set at 110 °C for a duration of 10 min.The obtained products underwent a triple washing procedure after natural cooling involving three cycles of rinsing with DMF and methanol, utilizing centrifugation (8000 r min −1 , 2 min) for each step.Finally, the specimens were desiccated in a vacuum oven set at 60 °C for the entire night.The synthesized samples were designated as 88C-1 in Table 1.Simultaneously, the synthesis conditions for 88C-2, 88C-3, and 88C-4 are documented in Table 1.Notably, the quantity of the organic ligand remained constant.

Fabrication of Carbon Coating Iron-Based Composites
The MIL-88C (Fe) prepared above served as the precursor and was then pyrolyzed under nitrogen at 800 °C for 2 h with a temperature rise rate of 5 °C min −1 .The obtained specimens derived from 88C-1, 88C-2, 88C-3, and 88C-4 were named as MD 1 , MD 2 , MD 3 , and MD 4 , respectively.

PPreparation of MOF-Derivatives/PVDF Coaxial Rings
The total mass of MOF-derivatives (MD 1 , MD 2 , MD 3 , or MD 4 ) and PVDF was 120 mg.The filling amounts of MOF-derivatives were 5 wt% (6 mg), 10 wt% (12 mg), 15 wt% (18 mg), and 20 wt% (24 mg), respectively.Initially, a certain number of MOF-derivatives and PVDF were dispersed by ultrasound in 5 mL of DMF for 30 min.Then, the homogeneous dispersion was transferred to an oven to evaporate at 90 °C for 4 h to obtain a flexible black film.The flexible black film was put into a special mold and pressed at 220 °C and 5 MPa for 15 min, then naturally cooled to room temperature to obtain the coaxial ring (Φ out = 7.00 mm and Φ in = 3.04 mm).

Material Characterization
The tests of crystal structure and phase composition of materials were conducted utilizing a Rigaku Smartlab SE X-ray diffraction (XRD) instrument employing Cu Kα radiation with a scanning rate of 5° min −1 .The surface features and internal constitution of the specimens were investigated via scanning electron microscopy (SEM, ZEISS Sigma 300) and transmission electron microscopy (TEM, FEI TECNAI F20).The elemental composition and spatial distribution were characterized via energydispersive X-ray spectrometry (EDS, Ultim Max 100).Raman spectra were documented using a high-resolution confocal microscopy laser Raman instrument (WITec alpha300R) with a 532 nm laser source.The elemental electronic structure and valence state of the materials surface were revealed via X-ray photoelectron spectroscopy (XPS, Thermo Scientific K-Alpha) employing Al Kα radiation.The magnetic hysteresis loops were tested by vibrating sample magnetometer (VSM, Lake Shore 7404).Thermogravimetry (TG, Netzsch STA 449 F5 Jupiter) analysis was carried out under N 2 atmosphere with a heating rate of 10 °C min −1 .Fourier transform infrared (FT-IR) spectra were acquired using an IRAffinity-1 spectrometer with spectral grade KBr sheets.The EM parameters of the coaxial rings were measured with a coaxial approach using a vector network analyzer (VNA Agilent TE5071C).

Characterizations of the As-Synthesized Materials
The preparation process of the derivatives from MIL-88C (Fe) and the absorbers mixed with PVDF is depicted in Scheme 1.The MIL-88C (Fe) has a three-dimensional (3D) hexagonal structure, in which trimeric units of Fe 3 O are linked by NDC 2-to form a periodic 3D structure of [Fe 3 O(NDC) 3 (H 2 O) 2 (NO 3 )] n [29,30].Anisotropic growth leads to a hexagonal rod shape (Fig. 1a).Interestingly, the M/O, oil bath temperatures, and oil bath times were able to change AR.Specifically, increasing the M/O, lowering the oil bath temperature, or extending the oil bath time can raise the AR of MIL-88C (Fe).Iron ions (Fe 3+ ) play a significant role in directional growth, and the larger their amount, the longer the nanorods grow.The temperature of the oil bath affects the nucleation rate.The lower the temperature, the fewer the crystal nuclei, and accordingly, the longer the nanorods will be.The prolongation of reaction time leads to a more complete growth of the crystal along the axial direction.In addition, the sample surfaces became rough, and some small crystals were observed on their surfaces when the oil bath time was long enough.This phenomenon reveals that dissolution-recrystallization may occur on the material surfaces.The pyrolysis process consists of two stages: carbonization and graphitization.For pure organic materials, the graphitization process requires high temperatures above 1000 °C.
If metal components are used as a catalyst, the temperature required for graphitization can be reduced to below 1000 °C [31].Therefore, in this work, iron and metal compounds containing Fe, generated during pyrolysis, play catalytic roles in the graphitization process.After pyrolysis, the derivatives shrank in size due to the decomposition of organic ligands.And their surfaces became rough with a large number of nanoparticles of variable sizes evenly distributed on surfaces, but they still retained the basic shape of the precursors (Fig. 1b/1c).In the case of MD 3 , for example, the results of EDS indicated the content of O was higher than it should be in Fe 3 O 4 (Fig. S1).It demonstrates that there are other forms of O present. Figure S2 illustrates the FT-IR results of MD 3 , with peaks observed at 1634 and 1400 cm −1 corresponding to the symmetric and asymmetric stretching modes, respectively, of the carboxyl groups in 2,6-NDC [32].Furthermore, the fact that the derivative retained the basic shape of the precursor also indicated that the organic ligand had not completely collapsed.Both of the above results suggest that the excess O is present in organic ligands.TEM and HRTEM were used to further obtain the component and structural information of materials (Fig. 2).For MD 1 , MD 2 , MD 3 , and MD 4 , the TEM results all showed that the iron-based components were encapsulated in carbon.The iron-based phases of MD 1 and MD 2 exhibited small and dense structures (Fig. 2a 1 /b 1 ), whereas MD 3 and MD 4 displayed bulky and sparse characteristics (Fig. 2c  The XRD result (Fig. S3) indicated the successful synthesis of precursor MIL-88C (Fe).As shown in Fig. 3a 1 , Fe-containing phases in the derivatives included Fe (JCPDS: 06-0696), Fe 3 O 4 (JCPDS: 88-0315), or Fe 3 C (JCPDS: 34-0001).For the carbon component, MD 1 and MD 2 were dominated by GC (JCPDS: 34-0567) (Fig. 3a 1 ).The broad peaks of MD 3 and MD 4 at about 20~30° represented the (002) plane of amorphous carbon (AC) (Fig. 3a 2 ) [33].The specific phases of the derivatives were as follows: MD 1 and MD 4 (Fe/Fe 3 O 4 /Fe 3 C@C), MD 2 (Fe/ Fe 3 O 4 @C), and MD 3 (Fe 3 O 4 @C), in good agreement with HRTEM and SAED results.Notably, Fe phases are present in all the composites except MD 3 .This occurrence may be attributed to the high M/O conditions, where the O content of the organic ligand is not sufficient to react with all Fe 3+ ions to form Fe 3 O 4 .Consequently, the surplus Fe 3+ ions are reduced to Fe phases.Raman spectroscopy, an effective bulk characterization method, is used for analyzing the carbonaceous materials.The peaks at the 1350 and 1580 Raman shifts in Fig. 3b represent the D and G peaks of the carbonaceous materials.The D peak serves as an indicator of imperfections within the carbonaceous materials, here including sp 3 -C in amorphous carbon and flaws in the sp 2 -hybridized mesh.The G peak indicates perfect graphite consisting of sp 2 -hybridized carbon [34,35].The I D /I G intensity ratio is frequently employed to depict both the level of graphitization and the presence of defects within carbon materials.Figure 3b shows the Raman spectra of the derivatives.The I D /I G results for MD 1 , MD 2 , MD 3 , and MD 4 were 0.81, 0.88, 1.03, and 2.82, respectively, indicating a decreasing order in the degree of graphitization.The aforementioned findings are also in alignment with the carbon diffraction peaks in the XRD results.The degree of graphitization was related to the type and size of the Fecontaining phases used as catalysts.Specially, the MD 2 had iron-based components with the same small sizes as MD 1 but lacked the Fe 3 C phase, which made it have a slightly lower graphitization degree than MD 1 .For MD 3 , it only had Fe 3 O 4 with a larger size and therefore had an elevated value of I D /I G as well as a lower level of graphitization.Although MD 4 had the same iron-based components as MD 1 , the larger size of the iron-based phases left less room for the carbon component and therefore also had a lower degree of graphitization.Moreover, the broad 2D peaks near 2700 cm −1 represent the presence of several layers of graphite, and the narrow 2D peaks represent multiple layers of graphite [36].The conclusions drawn from the 2D peaks were in accordance with the HRTEM results.The above results indicate that the content of Fe 3+ , the oil bath temperature, and the oil bath time affect the phases and degree of graphitization of the derivatives.
The XPS results are presented in Fig. 3c.The survey spectra (Fig. 3c 1 ) showed that the composites contained C, O, and Fe elements.The intensity of the C peaks indicated that MD 1 and MD 2 had higher and similar carbon content, whereas MD 3 and MD 4 exhibited sequentially lower carbon content.The peaks at 284.3 eV represented sp 2 -C, which demonstrated the material surface to be mainly graphite (Fig. 3c 2 ).Meanwhile, the peaks representing carboxyl groups from the undecomposed organic ligand 2,6-DNC were also present.In the narrow spectra of the O element (Fig. 3c 3 ), the two peaks with lower binding energies were attributed to the lattice oxygen (529.9 eV) and vacancy oxygen (Vo) in Fe 3 O 4 (531.6 eV) [37].The two remaining peaks at 532.7 eV and 533.5 eV denoted the oxygen in the carboxyl group, further demonstrating the existence of the undecomposed organic ligand 2,6-DNC.As shown in Fig. 3c 4 , the Fe element mainly consisted of two valence states of " + 2" and " + 3," which indicated that the material surface was mainly Fe 3 O 4 , while the Fe and Fe 3 C were mainly present in the bulk phases.
The magnetic properties of the composites are reflected by the saturation magnetization (M s ) and coercivity (H c ) from the VSM results (Fig. 3d).In particular, MD 4 presented a ferromagnetic hysteresis loop with the maximum M s of 92.5 emu g −1 and H c of 331.6 Oe, which is due to the fact that it contains multifarious and bulky ferromagnetic components and low carbon content.MD 3 displayed the minimal M s of 39.9 emu g −1 and H c of 131.5 Oe.This is attributed to the reason that it has only a single magnetic component Fe  https://doi.org/10.1007/s40820-024-01435-z© The authors ferromagnetic compositions, MD 1 has a higher carbon content than MD 4 .Therefore, the M s and H c of MD 1 are considerably lower than those of MD 4 .
To determine the suitable temperature for obtaining carbon-coated iron-based composites, the TG curve of 88C-3 under N 2 conditions was examined (Fig. S4).The mass loss was 4.4 wt% in stage I, corresponding to the loss of guest H 2 O molecules.The mass loss of 7.6 wt% in stage II was attributed to the removal of coordinated H 2 O and DMF molecules.With the increase in temperature, the large mass loss of 26.7 wt% in stage III corresponded to the decomposition of organic ligands.The mass loss of 12.7 wt% in the last stage IV indicated complete decomposition of organic ligands [38], culminating at ~ 800 °C with a total weight loss of 51.4 wt%.Consequently, derivatives with adequate pyrolysis can be obtained at 800 °C.

MA Performance of Absorbers Derived From MIL.88C (Fe)
The values of RL and EAB are major indicators to assess the MA properties.RL values can be counted using the subsequent two equations [39]: where the normalized input impedance of the absorbers and free space is represented by Z in and Z 0 , respectively.ε r and μ r denote the relative complex permittivity (ε r = ε′ − jε″) and complex permeability (μ r = μ′ − jμ″).f represents the incident EM wave frequency, d denotes the thickness of the absorbers, and c signifies the speed of light propagation in vacuum.RL < -10 dB means the absorber possesses MA properties, and EAB designates the frequency region wherein RL < -10 dB.To fully exploit the MA performance of the absorbers, we chose PVDF instead of wax as the matrix.The effect of PVDF and wax as substrates on the MA properties was compared using MD 3 as an example (Fig. S5).As shown in Fig. S5a 1/2 , the ε′ and ε″ of MD 3 /PVDF are significantly higher than those of MD 3 /wax for the same filling amount (10 wt%), consistent with previous studies [40].The notable increase in the dielectric constant is attributed to dipole polarization resulting from the difference in electronegativity between F/H atoms and C atoms, as well as interfacial polarization
Notably, the graphitization degree and AR together influence the optimum filler loading, with the degree of graphitization being the main factor.The greater the level of graphitization, the more readily a conductive network.MD 1 had the highest degree of graphitization and a large AR, resulting in an optimal filler amount as low as 5 wt%.MD 2 was slightly less graphitized than MD 1 but had a larger AR than MD 1 , so the optimum filler loading could also be as low as 5 wt%.Compared to MD 1 , MD 3 had a low degree of graphitization and a small AR; thus, it had an optimal filler amount of 10 wt%.Although MD 4 had the largest AR, the degree of graphitization was the lowest, therefore the optimum filler loading was 20 wt%.
Generally speaking, two essential factors dictate the MA properties.The primary consideration is impedance matching coefficient (M z ).For samples with finite thickness, utilizing M z to evaluate impedance matching characteristics is more reasonable [43], and value near 1.0 indicates a well-matched impedance, with a large number of EM waves able to be introduced into materials.

3
An assessment of the impedance matching capability of absorbers is conducted within the range of 0.8 to 1.0 [34].Besides good impedance matching, MA performance is also influenced by the attenuation constant (α).The value of M z and α were determined by the subsequent equations [34,44]: where the Z in ′ is the real input impedance, where the parameters ε′, μ′, ε′′, and μ'' are real and imaginary parts of relative complex permittivity and permeability, respectively.As shown in Fig. 5a 1 /b 1 , MD 2 /PVDF (5 wt%) possessed a greater impedance matching area delineated with the white dotted lines and more Z values close to 1 than MD 1 /PVDF (5 wt%).The values of α for both are close to each other, indicating similar attenuation capacities (Fig. 5c 1 ).As a result, MD 2 /PVDF (5 wt%) had a stronger RL and a wider EAB.In Fig. 5a 2 /a 3 , MD 1 /PVDF (10 wt%) and MD 1 /PVDF (20 wt%) had poor impedance matching.This is due to the high graphitization degree, which makes MD 1 /PVDF easier to form conductive networks, resulting in higher dielectric constants within high filler loadings (Fig. S8).Notably, for MD 1 / PVDF within a filler loading of 20 wt%, the ε′′ is greater than ε′, which indicates the formation of a conductive network resembling a conductor at high filler levels, leading to skin effects [45].Due to the impedance mismatch, even though MD 1 /PVDF (10 wt%) and MD 1 /PVDF (20 wt%) had higher loss capacities (Fig. 5c 2 /c 3 ), they could not exhibit MA properties.In contrast, both MD 3 /PVDF (10 wt%) and MD 4 /PVDF (20 wt%) had excellent impedance matching (Fig. 5b 2 /b 3 ), combined with suitable dissipation capacities, and therefore exhibited excellent MA performance.EM parameters were analyzed in depth for the purpose of further exploring the MA mechanism of the three excellent absorbers.The real parts (ε′ and μ′) signify the storage capability of electrical and magnetic energy, while the imaginary parts (ε″ and μ″) indicate the ability to dissipate EM energy.The ε′ values in Fig. S9a 1 exhibit an overall downward trend from 13 to 7, 8, or 10 approximately with increasing frequency, exhibiting the dispersion effect.The ε″-f graphs of the three absorbers were observed in Fig. S9a 2 , and the curves all exhibited an overall decreasing trend from 7.6 to 3.6, 5.2 to 3.3, and 5.7 to 2.9, respectively.This represents the conductivity loss and small resonance peaks represent the polarization loss [46].The MD 4 /PVDF had the highest number of resonant peaks, the MD 2 /PVDF the second most, and the MD 3 /PVDF the least.The dielectric loss tangent (tanδ ε = ε′′/ε′) was computed to assess the ability of absorbers to lose electrical energy, and the corresponding curves showed the same trend as ε″ (Fig. S9a 3 ).Figure S9b 1/2 display the μ′-f and μ″-f graphs of the three absorbers, and the presence of resonance peaks is observed, indicating that the absorbers have magnetic loss capability to some extent.Some negative μ″ values of the MD 2 /PVDF may result from the induced alternating magnetic field generated by the motion of free electrons within the alternating electric field, leading to the radiation of magnetic energy from the conductive networks [47].The values of the magnetic loss tangent (tanδ μ = μ″/μ) are also computed to assess the ability to lose magnetic energy (Fig. S9b 3 ).By comparison, the value of tanδ ε is higher than tanδ μ , demonstrating that the contribution of dielectric loss to the MA performance of absorbers predominates over magnetic loss.
Types of microwave loss can be classified from both component and micro/macroscopic geometric structure perspectives.The component level includes dielectric loss caused by dielectric materials and magnetic loss from magnetic materials.Dielectric loss comprises conduction loss and polarization relaxation within the range of 2 ~ 18 GHz.The Debye relaxation theory elucidates the mechanism underlying dielectric loss, and it can be summarized the correlation for ε′ and ε′′ as [48]: where ε s stands for the static dielectric constant as the frequency approaches 0, while ε ∞ denotes the permittivity in the limit as frequency approaches infinity.The ε′-ε′′ curve of the three absorbers exhibited a shape accompanied by several Cole-Cole semicircles in the high frequency range and a straight line in the low frequency range in Fig. S9c 1~3 .This indicates that the dielectric loss included the contributions of polarization and conduction losses.The semicircles depict the process of relaxation polarization, encompassing dipole polarization, defect-induced polarization, and interface polarization.The dipole polarization was caused by the variance in electronegativity between F/H atoms and C atoms in PVDF and dangling bonds on the surface of Fe 3 O 4 [49].The defect-induced polarization was generated by carbon defects [50], grain boundaries [51], and Vo in Fe 3 O 4 [52].The interface polarization was formed by various heterogeneous interfaces with different electronegativity, including GC, AC [53], PVDF, and iron-based components.The MD 4 /PVDF had the largest number of semicircles due to the most iron components (Fe, Fe 3 O 4 , and Fe 3 C), and these iron components produced more interfacial polarization.Additionally, the number of semicircles observed was in accordance with the regularity of resonance peak occurrences in the ε′′-f curves.The straight lines demonstrate the conductive loss attributed to GC.Generally, the magnetic loss within the 2 ~ 18 GHz range includes eddy current loss, natural resonance, and exchange resonance.The analysis of magnetic loss type involves the calculation of the value of C 0 using the following equation [54]: In Fig. S10, the three absorbers all exhibited a clear fluctuation trend between 2 ~ 8 GHz and a constant trend within the band of 8 ~ 18 GHz.This phenomenon proves that the magnetic loss was caused by natural resonance (2 ~ 8 GHz) and eddy current loss (8 ~ 18 GHz).As MD 2 / PVDF and MD 4 /PVDF contained more magnetic components, both fluctuated more significantly between 2 ~ 8 GHz, indicating more significant magnetic resonance losses occurred.
Besides the component perspective, the types of microwave loss in terms of micro/macroscopic geometric structure, including dissipation of the reflected EM wave and multiple scattering, also play a key role.The dissipation of the reflected EM wave can be analyzed based on λ/4 theory via the subsequent equation [55]: where t m represents the matching thickness, and f m represents for the frequency corresponding to the optimal RL value under the matching thickness.λ and c are the corresponding microwave lengths and the velocity of light in vacuo.|ε r | and |μ r | denote the moduli of the relative complex permittivity and permeability.The λ/4 model leads to a different phase of π between reflection waves from the front and back of absorbers, which results in a reflection wave extinction effect and promotes MA performance [55].The coordinate points corresponding to the actual situation of the three absorbers were compliant to the f m − t m curve within the error range (Fig. S11) [43].It demonstrates that optimum RL versus matching thickness is in accordance with a specific rule and the optimal RL values are obtained at the most matched thicknesses.Meanwhile, the M z values at optimum thicknesses are extremely close to 1.0.Furthermore, randomly distributed nanorods tend to scatter EM waves multiple times, extending the propagation path and facilitating the conversion of EM energy into heat.

MA Mechanisms of Derivatives From MIL.88C (Fe)
The potential mechanisms of MA for the absorbers are illustrated in Fig. 6.The combination of moderate dielectric constants and certain permeabilities facilitated optimal impedance matching and obtained suitable attenuation constants.The random distribution of graphite-containing derivatives constituted a 3D conductive network that promoted electron migration and hopping, thereby exhibiting conductive loss properties.Defects such as carbon defects, vacancies, and grain boundaries in the composites induced and trapped electrons, resulting in polarization induced by defects.Dipole polarizations were generated by inherent dipoles in PVDF and dangling bonds on Fe 3 O 4 surfaces.More interfacial polarizations were also created among multiple phases.The magnetic components Fe, Fe 3 O 4 , and Fe 3 C converted EM energy into thermal energy through magnetic resonance and eddy currents.All the mechanisms discussed above were examined at the component level.Furthermore, micro/macroscopic geometric structural aspects such as multiple scattering and λ/4 theory were also significant in MA.

Comparison of MA Performance
Eventually, a specific reflection loss (SRL) value is introduced to properly assess the MA capability considering both filler loading and thickness.It can be calculated as follows [56]: We compared the SRL values of the MD 2 /PVDF and MD 3 /PVDF with those of recent advanced absorbers, as shown in Fig. 7.The specific performance data of the above absorbers are listed in Table S1.Surprisingly, our work exhibited minimal SRL values compared to other absorbers, indicating stronger reflection loss, thinner thickness, and lower filler loading.

EM Simulation Results of the SGHS
To further investigate the absorption results of the SGHS, simulations were conducted using HFSS in accordance with the finite element method and transmission line theory.As shown in Fig. 8a, the SGHS model was created with the pink border representing an air box, the bottom blue plate serving as an aluminum plate, and the SGHS in the middle section.Two sets of master-slave boundaries were established (Fig. 8b, c), with the Freud excitation boundary acting as the incident direction and the scanning port serving as the corresponding detection direction (Fig. 8d).During the EM simulation, a periodic structure was assembled through repetition of the honeycomb model.Simultaneously, the EM parameters of MD 3 /PVDF (10 wt%) were imported into the software.EM waves were incident vertically from directly The MA performance of SGHS is determined by its structural features.Therefore, this study focused on investigating the impact of the parameters (H, r, and R) on the MA properties of the SGHS.Thickness (H) has vital importance in honeycomb structures and significantly influences MA properties [57,58].As illustrated in Fig. S12a, both the RL and EAB of the SGHS increased notably with an increase in thickness.The minimum RL min of -59.0 dB was achieved when H = 30 mm, r = 7.5 mm, and R = 9 mm. Figure S12b depicts the impact of r on the MA performance of the SGHS.RL strengthened as r decreased within a certain range, but continued reduction in r led to RL weakening.This is attributed to the decrease in r causing an increase in the gradient of the honeycomb holes, which, to some extent, promotes impedance matching, enabling greater penetration of EM waves from the upper half [59,60].Simultaneously, the significant gradient of the symmetric lower half facilitates better multiple scattering of EM waves, resulting in excellent MA performance.However, excessively small r hinders the entry of EM waves into the SGHS.Meanwhile, EAB increased with decreasing r, continuously expanding the effective absorption of the SGHS at low frequencies.Notably, when r (8.5 mm) and R (9 mm) are nearly equal, the structure resembles the conventional through-hole honeycomb structure, exhibiting poorer MA performance.This indicates a significant advantage of the SGHS over the conventional honeycomb structure.The impact of R on the MA performance of the SGHS aligned with that of r (Fig. S12c).The maximum EAB of 14.6 GHz (3.4 ~ 18.0 GHz) occurred at parameters H = 30 mm, r = 7.5 mm, and R = 13 mm.The gradient of the honeycomb hole is fixed when the difference between R and r is certain.Under the condition of fixed honeycomb hole gradients, the RL curves of the structures had the same trend when R and r were varied simultaneously (Fig. S12d).This demonstrates the substantial impact of the honeycomb hole gradient on the MA properties of the SGHS.In summary, the parameters of the SGHS can be chosen based on the particular demands for MA performance.
To further analyze the MA mechanism of the SGHS, Fig. 8 presents the distribution of electric field and vector, magnetic field and vector, volume loss density, and Poynting vector at 10 GHz.In Fig. 8f 1 , in the upper half of the SGHS, the electric field strength increases closer to the honeycomb hole, suggesting that a significant quantity of microwaves penetrates the honeycomb hole.The decreasing gradient of the electric field in the lower half of the SGHS indicates a gradual attenuation of microwaves during transmission.Figure 8f 2 shows that the local current mainly travels along the wall of the SGHS (y-direction), concentrating primarily above and below the cross section where the honeycomb holes are situated.The geometric properties of the SGHS wall enable a longer path for local current conduction, which is advantageous [57].The magnetic field distribution was approximated to the electric field distribution, and the results similarly indicated that massive microwaves entered through the honeycomb holes in the upper half of the SGHS and were absorbed in the lower half (Fig. 8g 1 ).The direction of the magnetic moments mainly aligned with the honeycomb wall (x-direction, perpendicular to the direction of the local current), and likewise, they predominantly concentrated near the cross section where the honeycomb holes were positioned (Fig. 8g 2 ).The results of the volume loss density distribution highlighted that the loss occurs throughout the SGHS (Fig. 8h).The Poynting vector, being an energy flow density vector in the EM field, visualized the distribution of EM energy in the SGHS (Fig. 8i).The results effectively validate the initial design conjecture.The gradient structure

Conclusion
In summary, we systematically synthesized a range of MIL-88C (Fe) precursors with varying AR under diverse preparation conditions, subsequently acquiring the corresponding carbon-coated iron-based derivatives using identical pyrolysis conditions.It is concluded that alterations in the preparation conditions of the MIL-88C (Fe) precursor impact the phases, graphitization degree, and AR of the derivatives, influencing the optimum filling amount, impedance matching, and MA properties.The MD 2 /PVDF, MD 3 / PVDF, and MD 4 /PVDF absorbers demonstrated outstanding MA properties with optimal filler loadings below 20 wt%, some as low as 5 wt%.The MD 2 /PVDF (5 wt%) achieved a maximum EAB of 5.52 GHz (1.90 mm).The MD 3 /PVDF (10 wt%) exhibited a RL min value of -67.4 dB at 12.56 GHz (2.13 mm).Additionally, the absorbers prepared in this study displayed minimal SRL values compared to most existing absorbers, signifying stronger RL, thinner thickness, and lower filler content.Through HFSS EM simulation, the integration of the material benefits of MD 3 /PVDF (10 wt%) with the structural advantages of the SGHS resulted in an EAB of 14.6 GHz and a RL min of -59.0 dB.It is well suited for the demands of real engineering applications and will prove to be promising candidates in the MA field.

Fig. 2 1
Fig. 2 1 TEM images, 2 ~ 3 HRTEM images, 4 SAED patterns, and 5 element mapping for Fe, O, and C of a MD 1 , b MD 2 , c MD 3 , and d MD 4 3 O 4 and high carbon content.The values of M s and H c for MD 1 (46.4 emu g −1 , 161.2 Oe) and MD 2 (44.5 emu g −1 , 181.2 Oe) are almost equal.The M s of MD 2 is slightly lower than that of MD 1 due to the absence of Fe 3 C in MD 2 compared to MD 1 .Although MD 1 and MD 4 have the same variety of

Fig. 6
Fig. 6 Schematic illustration of MA mechanisms for the derivatives of MIL-88C (Fe)

Fig. 7
Fig. 7 Comparison of SRL of MD 2 /PVDF and MD 3 /PVDF with the current advanced absorbers

Fig. 8 a
Fig. 8 a SGHS model; b, c two sets of master-slave boundaries; d sources of excitation for the model; e parameters related to the SGHS; f 1 Simulated electric field and f 2 electric field vector distribution of the SGHS at 10 GHz; g 1 Simulated magnetic field and g 2 magnetic field vector distribution of the SGHS at 10 GHz; h Volume loss density distribution of the SGHS at 10 GHz; i Poynting vector distribution of the SGHS at 10 GHz